Infrared study of the phonon modes in PrMnOs and CaMnOa 
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The infrared (IR) reflectivity spectra of orthoriiombic manganese perovskites PrMnOs and 
CaMnOa are studied in the frequency range of optical phonon modes at temperatures varying 
from 300 to 4 K. The IR phonon spectra of these two materials are analyzed by a fitting procedure 
based on a Lorentz model, and assigned to definite vibrational modes of Pnma structures by com- 
parison with the results of lattice dynamical calculations. The calculations have been performed in 
the framework of a shell model using short range Born-Mayer-Buckingham and long range Coulomb 
potentials, whose parameters have been optimized in order that the calculated Raman and IR active 
phonon frequencies, and lattice parameters match with their experimental values. We find a close 
correspondence between the values of the IR phonon frequencies of PrMnOs and CaMnOs, which 
shows that the substitution of the Pr^^ ions with Ca^^ results in a reduction of the frequency 
of medium- and high-energy IR phonons, and an increase of the frequency of those of low-energy. 
Nevertheless, the experimentally obtained IR phonon amplitudes of the two materials appear to 
be unrelated. A comparative study of the vibrational patterns of these modes reveals that most 
of them correspond to complex atomic vibrations significantly different from PrMnOs to CaMnOs 
which cannot be assigned only to a given type of vibration (external, bending, or stretching modes). 
In particular, these results confirm that the structure of CaMnOs is quite far from the ideal (cubic) 
perovskite structure. 

PACS numbers: 78.30.-j, 63.20.-e, 75.47.Lx 



I. INTRODUCTION 

Among the Mn oxides with perovskite structure 
[T] (the so-called manganites) , Pri_a;Caj;Mn03 is well 
known as the unique showing insulator behaviour over 
the whole range of doping (x) at zero magnetic field. This 
is usually assigned to a strong lattice distorsion, namely 
the tilting of MnOg octahedra (or an average Mn-O-Mn 
angle smaller than 180°) induced by the small size of 
the Pr and Ca ions, and leading to the localization of 
the itinerant carriers in the e^-orbitals. The phase dia- 
gram of Pri_2;Caa;Mn03 exhibits a wide variety of ground 
states EJ [1 H [5]. In particular for 0.3 < x < 0.75, 
a Ci?-type charge ordered state — where the Mn'^"'" and 
Mn^"^ ions are alternately arranged within the x-z planes 
(in the Pnma coordinate axes orientation) with a con- 
comitant ordering of the e^-orbitals — is stabilized below 
Tco = 250-200 K m ig and then a CE-type antiferro- 
magnetic spin arrangement takes place at T/y = 170-180 
K [2]. An interesting point is that Tjv is smaller than 
Tco- This suggests the existence of a strong electron- 
phonon interaction arising from the Jahn- Teller (JT) dis- 
torsions, while the exchange coupling between localized 
spins = I in the i2g-orbital is lower. Note also that 
at low temperature when 0.3 < x < 0.5, the CE-type 
spin ordering turns into a canted antiferromagnetic state 
{TcA < 100 K), which could consist of a mixture of fer- 
romagnetic and antiferromagnetic clusters External 
forces — like magnetic fields [31 S H [HI H [10], elec- 
tric fields [m [lU [131 [H [13 [IS], hydrostatic pressure 
[Tfl [18], and photon exposure [191 [20] — can destabilize 
the C£^-type charge/orbital and spin ordered state into 
the ferromagnetic metallic phase that is absent in the 



zero magnetic field phase diagram of Pri_a;Caa;Mn03. It 
is generally accepted that the competition between ferro- 
magnetism and antiferromagnetism is an important key 
to the colossal magnetoresistance effect of the mangan- 
ites. Such a phase competition results from the com- 
plex interplay of charge, orbital, spin, and lattice degrees 
of freedom. The latter is expected to be quite active 
in Pri„2.Ca2,Mn03 which is a strongly orthorhombically 
distorted manganites. 

Infrared (IR) spectroscopy is a very effective tool to 
study systems with strong electron-lattice coupling. Oki- 
moto et al. [21] measured optical spectra of a single crys- 
tal of Pri_j.Ca:rMn03 {x = 0.4) varying temperature and 
magnetic field. In particular, an increase of the number 
of optical phonons was observed at the charge/orbital 
transition, but their origin was not clearly established. 
Similarly phonon anomalies were observed at low tem- 
perature in the half-doped coumpound Pro.5Cao.5Mn03 
by Ta Phuoc et al. [22] . Classifying the optical phonons 
of the end-member compounds PrMnOs and CaMnOs 
may be particularly useful for further studies on the evo- 
lution of these modes in Pri-jjCa^^MnOs, and for identi- 
fying the additional phonon modes which should emerge 
in the charge/orbital ordered state. 

The structural determination of PrMnOa was obtained 
by neutron diffraction [231 [Ml [H] , and by X-ray diffrac- 
tion {T3[. Below - 900 K, PrMnOs with the Pnma or- 
thorhombic structure presents a strong JT cooperative 
distorsion of MnOg octahedra related to the ordering of 
the Bg-orbitals [26 , and undergoes at Tn « 95 K an 
A-type antiferromagnetic spin ordering |27j . The struc- 
ture of CaMnOs is also described with the Pnma space 
group at temperatures ranging up to at least 800 K |28j . 
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However, available structural data pS", "30] indicate that 
CaMnOa has an orthorhombically distorted perovskite 
structure much less pronounced than PrMnOa. This is 
a consequence of the substitution of Mn'^+ {tg ej) with 
Mn''^"'" (tg) that actually induces two effects of different 
origins. One arises from the fact that Mn'*+ has no or- 
bital degeneracy, and therefore does not produce any JT 
distorsion. The other is due to the mismatch between the 
ionic radii ru of the i?-ions (Pr^"*" or Ca^"*") and the ionic 
radii rMn of the Mn'^+ or Mn^+ ions, which is measured 
by the tolerance factor t = {rR + ro)/V^{rMn + tq) 
(where t — 1 corresponds to an ideal perovskite struc- 
ture, whereas t <1 indicates an average Mn-O-Mn angle 
that gets smaller than 180°). Gd?^ and Pr'^^ having sim- 
ilar ionic radius and those of Mn*"*" (0.53 A) and Mn'^+ 
(0.645 A) being significantly different, the tolerance fac- 
tor of PrMnOa (0.94) turns out to be smaller than that 
of CaMnOa (~ 1), which should appear to have a struc- 
ture very close to the ideal cubic perovskite. Recently, 
the Raman spectra of PrMnOg ED [Ml 1331 131 ISS] 
and CaMnOg [35] Hi [371 [H have been studied in 
detail. In particular, it appeared that the two dominant 
peaks in the Raman spectrum of PrMnOa are activated 
by stretching of the basal oxygen ions (02) of the MnOe 
octahedra due to the JT effect, in striking contrast with 
the spectrum of CaMnOs where the two main peaks are 
activated by rotation of the octahedra around the [101] 
cubic axis. Similarly, IR-active phonon modes in these 
materials are expected to be very sensitive to different 
types of distorsion (rotations of MnOg octahedra around 
the [101] or [010] cubic directions, JT distorsion, and 
shift of R atoms) [35^. However, there has not been any 
systematic study of IR-active phonons in PrMnOa, at 
least to our knowledge. Concerning CaMnOa, reflectiv- 
ity spectra of Lai_3,Ca2.Mn03 have been measured from 
5 meV up to 30 eV for different values of x including 
X — 1, but only the electronic part of these spectra was 
discussed |40]. Furthermore, Fedorov et al. [41] have 
presented IR absorption spectra of CaMnOa, but the as- 
signement of the phonon peaks was done by comparison 
with the IR absorption spectra of LaMnOa which is a JT 
compound. 

In this paper, we present an analysis of IR-active 
phonons in PrMnOa, and CaMnOa. The experimentally 
observed phonon peaks are assigned to definite atomic 
vibrations by comparison with the results of lattice dy- 
namical calculations (LDC) for both these materials. 



II. EXPERIMENTAL DETAILS 

Very fine powders of high chemical homogeneity of 
PrMnOa and CaMnOa were prepared using a modified 
citrate process as described in Ref.SH The polycrys- 
talline sample of PrMnOa was obtained from these par- 
ticles of PrMnOa pressed isostatically under 3 kbar in 
the form of a cylindrical slab of about 1 cm in diame- 
ter and 3 mm in thickness, and sintered at 1250°C in 



air for 12 hours. The polycrystalline sample of CaMnOa 
was synthesized in the same way, but it was sintered at 
1300°C in air for 24 hours, and then annealed at 800°C in 
oxygen current flow for 1 hour. Powder X-ray diffraction 
profiles recorded at room-temperature revealed that both 
samples had single-phase orthorhombic Pnma structure 
with the lattice parameters a = 5.7852 A, b — 7.5939 A, 
c = 5.4490 A for PrMnOa, and a = 5.2785 A, b = 7.4435 
A, c = 5.2585 A for CaMnOa. These data are in good 
agreement with the lattice parameters of stoichiometric 
samples reported previously p3 | Ei ) [25] [5111 H3] . 

Another indication of the quality of our samples was 
obtained from Raman measurements taken at room- 
temperature, which showed that both samples have Ra- 
man spectra very similar to those reported in the litera- 
ture 1311 IMl 1331 El 1351 131 133 IMllSn]. 

We have measured near normal incidence reflectivity 
spectra with a Bruker IFS 66v/S in the range 50 cm"^- 
17000 cm"-'^ from 4 K to 300 K with a continuous flow 
cryostat. The samples were polished up to an optical 
grade of 0.25 /Ltm. After the measurement, the samples 
were coated with a silver film and used as reference to 
calculate the absolute reflectivity. 

The real part of the optical conductivities (Ji{uj) were 
obtained from the reflectivity data by means of Kramers- 
Kronig transformations. In order to perform this calcula- 
tion, the reflectivity was simulated by a Lorentz oscillator 
centered around 24000 cm^^ in the high-frequency range, 
and set at a constant value between and 50 cm~^. We 
verified that such an assumption produces no noticeable 
effect on the conductivity spectra presented below. 

LDC were carried out by means of a shell model using 
the general utihty lattice program fGULP)|ii]. 

III. RESULTS AND DISCUSSION 
A. Spectra of PrMnOa and CaMnOs 

Figs.jlja) and 1(b) respectively display the real part of 
the optical conductivity of PrMnOa and CaMnOa at 300 
and 6 K in the far-IR frequency range which corresponds 
to the region of optical phonon vibrations. These spectra 
are consistent with the optical conductivity typical of an 
insulator that shows several peaks associated with the 
IR-active transverse-optical modes. As the temperature 
is lowered, no appreciable spectral change is observed, 
and the phonon features become narrower and sharper, 
as it is usually expected. 

The spectra shown in Figs, [ija) and 1(b) appear to 
be rather qualitatively similar in structure. They ex- 
hibit, around 160, 300, and 550 cm^^, three main phonon 
groups showing a fine structure. In particular, it may 
be seen that the phonon group located between 250 
cm~^ and 375 cm~^ — known as the bending band — 
is split for both PrMnOa and CaMnOa, in agreement 
with what is expected in orthorhombically distorted per- 
ovskites with Pnma structure |45j . The fine structure in 
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FIG. 1: (Color online) Infrared phonon modes in (a) PrMnOs 
and (b) CaMnOa, at 6 K and 300 K. The red solid lines 
represent the best fit of Eq.[l] to the low-temperature data. 
For better clarity, the conductivity values have been shifted 
by (a) 75 f2~^cm~^ and (b) by 150 f2~^cm~^, for the room 
temperature curves. 



this band is expected to be very sensitive to basic dis- 
torsions -D[ioi] and -D[oio] associated respectively to tilts 
of the MnOe octahedra around the [101] and [010] cubic 
axes. So, since Z)[ioi] (£'[010]) i^i PrMnOa is 1.26 (1.34) 
times larger than in CaMnOa |3TJ|38], the fine structure 
of the bending band in CaMnOa should appear signif- 
icantly reduced compared to that of PrMnOa. Indeed, 
such a general trend is clearly seen in Figs.[l];a) and 1(b). 

On the other hand, the high-energy modes (centered 
around ^ 550 cm"-"^) are thought to involve mainly 
stretching vibrations of MnOg octahedra. The frequen- 
cies of these modes are therefore expected to be directely 
related to the force constants of Mn-0 bonds, which are 
determined by the interatomic distances. Taking into 
account the much shorter Mn-0 distances in CaMnOa 
compared to PrMnOs [151 UHl 135] , one could expect 
the peak positions of the stretching modes to satisfy the 
inequality uocmo > ^pmo- As seen from Figs. [T| and |2] 
this is obviously not the case. It is worthy of note that 
the high-energy IR modes of LaMnOs — in which the 
Mn-0 distances are also greater than those in CaMnOa 
— were also observed at higher frequencies than those of 
CaMnOa in a previous study of the IR absorption spectra 
^41j . Such a behaviour can be explained by the fact that 



the charge replacement of Mn^+ in PrMnOs (LaMnOa) 
with Mn^"*" in CaMnOs is accountable for the decrease of 
the frequencies of these high-energy IR modes. 

The low-energy phonon group — corresponding to ex- 
ternal modes — mainly involves i? = Pr or Ca atoms. 
Therefore, the frequencies of these modes depend on the 
mass Affl of the atom i?, and are sensitive to force con- 
stants of R-0 bonds, which are determined by the i?-0 
distances dp-o- In an harmonic oscillator approximation 
{lo = y/k/Mji), assuming a spring constant k governed 
by l/{dH-o)^ dependence [46 , and using the values of 
mean R-0 distances available in Refs.[25]and[35l the ex- 
ternal mode frequencies in CaMnOa and PrMnOa should 
scale as ojcmo/^pmo — 1-9. Figs. [l] and |2] obviously 
show that this is not in quantitative agreement with the 
experimental data. So, such a simple interpretation in 
terms of harmonic oscillator is inadequate to explain the 
frequency shift of these modes, and one has to take into 
account the more complex orthorhombic phonon vibra- 
tions. Finally, a comparison of Figs.jlja) and 1(b) reveals 
very strong attenuation of the amplitude of some of these 
modes in CaMnOa. 

It is worthy of note that the present data which show 
a high number of observed phonons for CaMnOa and 
some similarity above ~ 220 cm~^ between the spectra 
of PrMnOa and CaMnOa, agree with the above men- 
tioned study of the IR absorption spectra of LaMnOa 
and CaMnOa by Fedorov et al. [IT]. The existence 
of such a correlation between the spectra of these per- 
ovskites shows that the Djt basic lattice distorsion — 
which originates from the JT effect in the case of PrMnOa 
and LaMnOa — is not the driving force of the orthorhom- 
bic distorsion in both these materials since their respec- 
tive Djt distorsion is almost 60 times larger than in 
CaMnOa [SI [35] • 



B. Oscillators model Fit 



To examine in some detail the phonon modes, we have 
fitted the spectra shown in Figs.[lja) and 1(b) using the 
well known Lorentz model: 



^)2+^2^2 



(1) 



where Sj, jj, and luj are respectively the amplitude, 
width, and frequency of the jth oscillator. The red solid 
lines in Figs.[l]display the best fit obtained using the min- 
imum number of oscillators that allow to achieve a reli- 
able fit to the spectra at low temperature. It is convenient 
to define the normalized amplitude Sj = 100 x Sj / Sj 
that corresponds to the percentage of spectral weight as- 
sociated with the jth oscillator. We list in Table |T] the 
values of the fitting parameters. 
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TABLE I: Best-fit values of tlie normalized amplitudes (Sj), widths (7j), and frequencies {coj) of the low temperature spectra 
shown in Figs.[T][a) and 1(b). The frequency and the symetry of the different phonon modes obtained by means of the lattice 
dynamics calculations (LDC) are also reported, ujj and ■jj are in cm~^, Sj is in %. 



PrMnOs CaMnOa 





Best fit 

7i 




LDC 

Symmetry 






Best fit 

Ij 




LDC 

Symmetry 




1.6 


21 


115 




125 








B2u 


133 


0.7 


17 


141 




132 


0.6 


17 


166 


Bzu 


172 


2.2 


22 


162 


B2u 


148 


1.1 


8 


176 


Bsu 


181 


6.1 


13 


175 


Biu 


184 


0.7 


8 


187 


Biu 


190 


6.9 


6 


184 


Bsu 


197 


0.8 


8 


227 


Biu 


224 


7.6 


17 


244 


B^u 


241 








B2u 


248.5 


9.2 


26 


257 


Biu 


279 


5.6 


16 


254 


Bsu 


250 


7.6 


28 


283 


Bsu 


298 


20.0 


18 


276 


B2u 


275 


1.2 


24 


298 


Bsu 


301 








Bzu 


281 


20.9 


24 


333 


Bxu 


348 


16.5 


20 


266 


Biu 


284 


7.1 


18 


349 


B2u 


351.5 


2.2 


12 


286 


Biu 


298 


0.9 


16 


369 


Biu 


370 


3.9 


16 


317 


B2u 


324 


3.8 


16 


360 


B2u 


379 


8.7 


13 


340 


Bsu 


345 








Bsu 


387 


6.4 


14 


332 


Biu 


346.5 


0.6 


10 


412 


Bsu 


396.5 


1.1 


13 


361 


Bsu 


364 


1.1 


17 


423 


Biu 


403 


2.2 


12 


348 


B2u 


368 


1.2 


18 


433.5 


B2u 


421 








Biu 


379.5 


1.1 


12 


444 


Biu 


426 


3.3 


12 


412 


Biu 


414 


1.3 


22 


480 


Bsu 


486 


0.2 


14 


449 


B2u 


450 


1.8 


14 


515 


Biu 


517.5 


1.2 


17 


457.5 


Bzu 


460 


0.9 


14 


524 


B2u 


525 


2.2 


14 


494 


Biu 


496 


2.8 


17 


557 


Bsu 


555 


2.1 


15 


510 


B2u 


504 


10.5 


18 


558 


B2u 


559 


9.2 


15 


539 


Bsu 


534 


2.4 


16 


571 


Biu 


562 


10.5 


16 


548 


Biu 


544 


0.5 


10 


580 


Bsu 


584 


1.5 


36 


584 


Bsu 


613 



C. Lattice dynamical calculations 

In order to identify the observed phonon peaks, LDC 
for both PrMnOa and CaMnOa have been performed by 
means of a shell model using GULP [H] . In this model, 
each ion is represented by a massless shell of charge Y 
and a core of charge X which are coupled by an har- 
monic spring constant K. The interionic interactions are 
described by long range Coulomb potentials and short 
range potentials of the Born-Mayer-Buckingham form: 

V{r) = A exp{-r/p) - C/r^ 

where r is the interatomic distance. These two-body 
radial potentials are of course a first approximation in 
which the effect of non-spherical potential terms of the 
interactions and other corrections have been neglected. 
However, it was shown that LDC using a shell model with 
this type of interionic interaction potential yield valuable 
information about the phonon eigenmodes of various oxy- 
des and can be used to perform the assignment of their 
experimental phonon spectra |47j . 

Starting from experimental crystallographic data, the 
lattice energy is minimized with respect to unit cell vari- 
ables — namely, the lattice parameters and atomic co- 
ordinates — while large potential parameter space zones 



are explored. We determine the potential parameters 
such that the lattice energy is minimum, and they repro- 
duce the phonon frequencies close to those determined 
by Raman and IR spectroscopy experiments. The values 
of parameters A and p for the 0-0 short range inter- 
action are taken from Refs.HS] and US] which have been 
successfully used for the modeling of oxides [371 [SUj- For 
PrMnOs and CaMnOa, we obtain the best fit with the 
parameters listed in Table [iTl It is worthy of note that the 
values of the potential parameters for the Mn'^+ ions de- 
rived in the present work are very close to those obtained 
by Islam et al. [5T] . 

To check the validity of our LDC, we show in Ta- 



ble |IIT] a comparison between the calculated structural 
parameters and those observed experimentally |23l |24l 
US [Ml US ISO]- It can be inferred that the LDC per- 
formed here produce a reasonable agreement with the 
experimental crystallographic data. In particular, the 
calculated lattice constants are very close to the ob- 
served values with a discrepancy of less than 2%. On 
the other hand, in accordance with the results of the 
group theory for the F-point of the orthorhombic Pnma 
structure of i^MnOs [IH US [53], we obtain 60 vibra- 
tion modes of which 24 {7 Ag + 5 Big + 7 i?2g + 5 B^g) 
are Raman active, 25 (9 -I- 7 B2u + 9 B^u) are IR- 
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TABLE 11: Sell model parameters. Cut-off radii of 10 A and 12 A were respectively used for the cation-oxygon and oxygen- 
oxygen short range potentials. 



Ion 


X{\e\) 


Y(\e\) 


K (eV A~2) 


Ionic pair 


A (eV) 




C (eV k-"") 


Pr 


1.5 


1.5 


30 


Pr-0 


45461 


0.2272 





Ca 


-0.6 


2.6 


80 


Ca-0 


10040 


0.2401 





Mn^+ 


0.5 


2.5 


95 


Mn3+-0 


1337 


0.3260 





Mn*+ 


1.9 


2.1 


40 


Mn*+-0 


4980 


0.2650 





O 


0.24 


-2.24 


42 


O-O 


22764 


0.1490 


27.88 



active, 8 (8 A^) are silent, and 3 {Bi^ -t- + B^^) are 
acoustic. For both PrMnOs and CaMnOa, we find good 
agreement between the calculated Raman-active phonon 
modes and the available experimental Raman spectra 
[nilSUIMllSllSSllMllSZllMllSn]- in particular, the av- 
erage deviation between the calculated and observed Ra- 
man frequencies is less than 10%. Indeed, the two domi- 
nant peaks observed by Raman spectroscopy in PrMnOa 
around 605 and 485 cm~^ and assigned to the symmetric 
[i?2g(l)] and asymmetric [j4g(I)] stretching of 02 ions in 
the xz planes [3TJ [32l |33l [34l [35] are respectively pre- 
dicted by the present calculation at 610 and 466 cm~^ 
(the numbering of the modes follows Refs.[2IJ ESI I3S1 1551 
and [5^. Moreover, the other main Raman peaks ob- 
served in PrMnOs around 465, 450 and 325 cm~^ which 
have been respectively assigned |3I1 EH IMl Ell ES] to the 
MnOe bending [Ag{i) + B2g{i)\ and tilt [Ag{4)] modes, 
are reproduced by our calculation at frequencies 510, 508 
and 357 cm~^, respectively. In the CaMnOa case, it is 
admitted [351 [SH] that the two pairs of Raman peaks ob- 
served around 487 and 466 cm^^ and around 258 and 
245 cm~^ cannot be related to Ag(I)+_B2g(I) (stretching) 
modes — the dominant ones in PrMnOa. The Raman fre- 
quencies of the doublet around 487 and 466 cm~^ which 
has been assigned [33 [35] to Ag{i) + -629 (3) (bending) 
modes, were obtained from a previous shell model com- 
putation by Abrashev et al. [38] at 467 and 453 cm~^, 
whereas the present calculation gives these frequencies at 
510 and 503 cm~^. Finally, it should be emphasized that 
the assignment of the pair around 258 and 245 cm~^ still 
remains controversial. By comparison between observed 
and calculated Raman frequencies f35^, this doublet was 
related to Ag{7) + -629 (7) modes involving the displace- 
ments of the apical oxygen (01) and Ca atoms, while 
— based on the correlation between frequencies and rel- 
evant structural parameters of Lai-j^CajjMnOs series — 
Martm-Carron et al. [3S] assigned it to Ag{A) + _B2g(4) 
(tilt) modes. 



D. Assignment of infrared-active phonon modes 

As our measurements were done on polycrystalline 
samples, the IR-active phonon peaks corresponding to 
modes of Bin, B2U and symmetry are undistinguish- 
able. For this reason, we make the assignment in order 
that the average deviation between the calculated and fit- 




FIG. 2: (Color online) Correspondence between the fitted 
(calculated) infrared active-phonon frequencies of PrMnOa 
[abscissae axis] and those of CaMnOs [ordinates axis]. The 
blue (red) symbols denote the fitted (calculated) frequencies; 
the meaning of filled symbols is given in the text, and the 
straight line is a guide for the eyes. Since the infrared spectra 
are obtained from polycrystalline samples, the assignment of 
experimental peaks to Bin, B2U and B^u modes is made by 
comparing the calculated frequencies with those determined 
from the best fit of the Eq.[l]to the data. 



ted phonon frequencies is minimum. The results of this 
comparison are shown in Table [l] It can be noted that 
very few ambiguities about the assignment of phonons re- 
main, such as, for instance, the two peaks in the spectrum 
of PrMnOa at 360 and 369 cm~^ which could be equally 
related to the Biu and B^u modes at 370 and 379 cm~^, 
respectively. In such cases, the assignment of the peaks 
is obtained by comparing the experimental amplitudes to 
those estimated qualitatively within the framework of a 
rigid-ion approach, as described below. Taking into ac- 
count the assignment of the modes and representing the 
fitted (calculated) values of ojj by blue (red) symbols, we 
show in Fig.|2] the relationship between the IR phonon 
frequencies of PrMnOa and CaMnOa . Note that the ex- 
perimental frequencies at 162, 244, 283, and 423 cm^^ 
for PrMnOa do not appear since they have no counter- 
part in the fitted spectrum of CaMnOa, and vice versa 
for the modes at 340 cm~^ for CaMnOa (in all these 
cases, the calculated frequencies are represented by filled 
symbols in the figure). It is clear from Fig. |2] that the IR 
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TABLE III: Comparison of the lattice parameters of PrMnOa and CaMnOs observed experimentally at room temperature 
{Pnma space group) with those obtained by means of the lattice dynamics calculations (LDC). The last seven lines represent 
the fractional atomic coordinates; R (Pr/Ca), Mn, apical (01) and basal (02) oxygen atoms are respectively located at Wyckoff 
positions 4c {x, 1/4, z), 4b (0, 0, 1/2), 4c {x, 1/4, z) and 8d {x, y, z), and at their symmetry equivalent positions. 

PrMnOs CaMnOa 





Observed" 


Observed" 


Observed" 


Observed' 


LDC 


Observed' 


Observed' 


Observed' 


LDC 




Ref. EH 


Ref.El 


Ref.Eg 


This work 




Ref.El] 


Ref.[29] 


Ref.;30li 




a (A) 


5.5914 


5.786 


5.8129 


5.7852 


5.6658 


5.28287 


5.2813 


5.279 


5.3110 


6(A) 


7.672 


7.589 


7.5856 


7.5939 


7.7042 


7.45790 


7.4582 


7.460 


7.4454 


c(A) 

V (A3) 


5.4562 


5.450 


5.4491 


5.4490 


5.4739 


5.26746 


5.2687 


5.273 


5.2446 


234.05 


239.31 


240.278 


239.386 


238.944 


207.5332 


207.528 


207.657 


207.388 


R (4c) X 


0.0451 


0.066 


0.0639 


0.062 


0.0517 


0.0330 


0.0313 


0.030 


0.0320 


z 


0.9931 


0.983 


0.9911 


0.991 


0.9913 


0.9944 


0.9929 


0.995 


0.9924 


Ol (4c) X 


0.4825 


0.484 


0.4819 


0.482 


0.4708 


0.4990 


0.4869 


0.491 


0.4701 


z 


0.0770 


0.084 


0.0834 


0.105 


0.0757 


0.0657 


0.0624 


0.066 


0.0782 


02 (8d) X 


0.2978 


0.313 


0.3174 


0.311 


0.2920 


0.2870 


0.2882 


0.288 


0.2876 


y 


0.0404 


0.042 


0.0430 


0.042 


0.0404 


0.0332 


0.0342 


0.033 


0.0417 


z 


0.7156 


0.716 


0.7151 


0.730 


0.7058 


0.7116 


0.7118 


0.709 


0.7112 



"by neutron diffraction 
'by X-ray diffraction 



phonon frequencies observed in PrMnOs and CaMnOa 
are closely correlated, and that this result is quantita- 
tively in good agreement with the present LDC. On the 
other hand, we also made a comparison between the fit- 
ted values of Sj for each IR phonon mode of PrMnOa 
and CaMnOa , and no correlation between them has been 
found. This suggests that the replacement of the Pr'^+ 
ions with Ca^"*" drastically changes the vibrational pat- 
terns of IR phonons. 

E. Analysis of the infrared active modes 

Let us briefly remember some general considerations 
about the expected forms of the IR-active modes [54] for 
discussing — on the basis of the present assignment — 
the atomic vibrations related to the phonon peaks ob- 
served in PrMnOa and CaMnOa. Since the Pr/Ca and 

1 atoms in the Pnma space group occupy Wyckoff posi- 
tion 4c, their motions are restricted by the site symmetry 
(Cg^) within the xz planes for the and B^u modes, 
and along the y direction for the mode. As for the 
Mn and 02 atoms — which are respectively localized at 
Wyckoff positions 46 (or 4a) and 8d — there are no sym- 
metry restrictions on the direction of their vibrations for 
all IR-active modes. Nevertheless, we can show using the 
results of the group theory that the general movement of 

02 atoms can be reduced for each mode of symmetry 
Biu, B2u, and B^u to a combination of three indepen- 
dent displacements involving in-phase or out-of-phase vi- 
brations along (stretching) or perpendicular (bending) 
to the Mn-02 bonds. The shape of these 02 vibrations 
is shown in Fig.|3] We see that the B2u and i?3„ 
bending vibrations of 02 atoms can be either in-phase 
or out-of-phase, while the Bi^ and B^^ stretching vi- 
brations are in-phase and those with symmetry B2u are 



in-phase out-of-phase 

stretching 

bending bending 




FIG. 3: (Color online) Classification of vibrations of basal 
oxygen (02) atoms for each mode of symmetry B\u, B2u, 
and B^u- Because the unit cell of RMnOa compounds with 
the Pnma structure contains two MnOe octahedra along the 
y direction, the 02 atoms in these neighboring octahedra may 
vibrate in-phase or out-of-phase. 



out-of-phase. Moreover, it can be noted that the _Bi„ and 
B^u in-phase bending vibrations of 02 atoms are related 
to one another by a rotation of 90° of MnOe octahedra 
around the y axis. The same occurs with the Biu and 
B^u stretching vibrations of 02 atoms. 

Taking into account the effects of the atomic polariza- 
tion, an expression for Sj is given in Refs.lM] and [Ml In 
particular, for every orthorhombic crystal one finds that 
Sj = |Pa(i)| where f2 is the volume of the unit cell, 
and Paij) is the nonzero component of the electric dipole 
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FIG. 4: (Color online) Atomic displacements of low-energy modes of Biu, B2u, and Bsu symmetry of PrMnOs (top) and 
CaMnOa (bottom) with Pnma structure. The numbers in brackets correspond to the calculated frequencies in cm~^. 



moment vector of the jth IR-active mode, namely Pz for 
Biu, Py for B2u, and p^ for B^u in the Pnma coordi- 
nate axes orientation. If the shell vibrational contribu- 
tion is neglected, then it can be shown [55j |56j |57] that 
PaU) = X]fc "^kU) where the subscript k labels the 
atoms in the unit cell, eZk is the effective charge of the 
Kth atom, and u"{j) denotes its active displacement com- 
ponent in the jth mode, so that Sj = IjEk u"{j)]'^. 
For instance, as can be seen in Fig.|3] the out-of-phase 
bending vibrations of 02 atoms (for all symmetries) and 
the stretching vibrations of 02 atoms with i?2u symmetry 
produce a zero dipole moment. Of course, such a formula 
for Sj derived from the rigid-ion model is very rough. 
However, for our purposes an important feature of this 
result is to highlight the possibility that some IR-active 
modes characterized by relatively important atomic vi- 
brations have a small resultant dipole moment, and hence 
a small oscillator strength Sj. As seen in Fig.ffl the 
low-energy modes of CaMnOa are the examples of such 
modes. Indeed, each of these modes does not produce 
a large dipole moment. For instance, the z component 
of atomic displacements giving the dipole moment in the 
Bill modes is small for the one at 187 (190) cm^^ (in 
the following, the numbers in brackets indicate the cal- 
culated frequencies) and is not represented in Fig. |4] The 
same occurs with the x component for the B^u modes at 
166 (172) and 176 (181) cm^^. The Bi„ and B2„ modes 
at 227 (224) and (133) cm^^ provide another example of 
modes carrying a small dipole moment. For these modes, 
the active components of the Ca and Mn vibrations are 
relatively large but in the opposite direction while those 
of the 01 and 02 atoms are small, which leads to a small 



resultant dipole moment. In PrMnOs, the situation is 
clearly different with the low-energy IR-active modes. At 
first glance, the Bi„ and modes at 175 (184) and 162 
(148) cm~^ are respectively similar to those obtained at 
227 (224) and (133) cm^^ in CaMnOg but the active 
components of the 01 and 02 vibrations are large and in 
the same direction, and hence, the dipole moment of both 
these modes is not small. Furthermore, we also find that 
the dipole moment of the B^^ mode at 184 (197) cm~^ is 
not small due to a large x component of atomic displace- 
ments. So, our calculations predict that unlike PrMnOa, 
none of the low-energy IR-active modes of CaMnOa has 
a large oscillator strength. It can be clearly seen from 
Fig.[T] that this prediction is in reasonably good agree- 
ment with the experimental data. 

Fig.[5]shows a representative selection of phonon modes 
of PrMnOa involving almost no movement of Pr atoms, 
as expected for internal modes. As can be seen, the 
modes at 333 (348), 480 (486), and 558 (559) cm"! — as 
well as those at 283 (298), (387), 412 (396.5), 423 (403), 
and 515 (517.5) cm~^ — correspond to complex motions 
of MnOe octahedra for which the amplitudes of the vibra- 
tions of Mn atoms are comparable with those of oxygen 
atoms. For the mode at 580 (584) cm~^ (see Fig.^ — 
as well as those at 360 (379), 433.5 (421), and 444 (426) 
cm~^ — only the 02 and Mn atoms play a role, so that 
the 01 atoms can be assumed to be at rest. On the 
other hand, there are also modes where only the oxygen 
atoms vibrate. In the mode at 571 (562) cm~^ displayed 
in Fig. [5] — as well as the 557 (555) cm~^ mode — the 
01 and 02 atoms vibrate with comparable amplitudes, 
whereas for the mode at 524 (525) cm~^ (see Fig.^ the 
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vibrations of 02 atoms dominate. As was stated above 
and shown in Fig. [3] the general movement of 02 atoms 
is a combination of bending and stretching vibrations. 
The Biu mode at 571 (562) cm~^ provides an example 
of such a mixture consisting of ^80% of stretching vibra- 
tions, ~17% of in-phase bending vibrations, and ^3% of 
out-of-phase bending vibrations. Similarly, the mode 
at 524 (525) cm^^ consists of ^^88% of stretching vibra- 
tions, ^--^11% of out-of-phase bending vibrations, and less 
than 1% of in-phase bending vibrations. As precisely this 
latter type of vibration produces the dipole moment, its 
small magnitude leads to a very weak intensity. It should 
be emphasized that this mode is also found by Smirnova 
[53] in its LDC of LaMnOs (-82^ mode at the frequency 
of 562 cm~^). On the other hand, there are also modes 
for which the movement of 02 atoms correspond to ei- 
ther quasi-pure bending vibrations or quasi-pure stretch- 
ing vibrations. So, the modes at 480 (486), 558 (559), 
and 580 (584) cm~^ can be classified as pure stretch- 
ing, as pure in-phase bending, and as pure out-of-phase 
bending, respectively. It is interesting to remark that the 
in-phase bending mode at 558 (559) cm~^ — which 
is characterized by a large dipole moment resulting from 
displacements of oxygen and Mn atoms along the y axis 
— is also obtained in the LDC of LaMnOs reported by 
Smirnova [53] {B2U mode at the frequency of 516 cm~^). 

Another significant point is that the atomic displace- 
ment patterns of the mode at 571 (562) cm"-'^ (see 
Fig. [5]) is linked to that of the B3U modes at 557 (555) 
cm by a rotation of 90° around the y axis. The same 
occurs for the -Bi„ and B^u modes at 423 (403) and 412 
(396.5) cm"-'^. For both of these pairs, the frequency dif- 
ference between the corresponding _Bi„ and i?3„ modes 
is expected to be small, and very sensitive to orthorhom- 
bic distortions. The experimental values A/ « 11 — 14 
cm~^ turn out to be large compared with the calculated 
values of 7 cm^^. This suggests that the present LDC 
does not fully account for the strutural distortion effects 
in PrMnOa. 

Fig. [6] shows a representative selection of the phonon 
modes of CaMnOa. As can be seen, the modes at 266 
(284), 340 (345), and 361 (364) cm-\ as well as those at 
(248.5), 254 (250), (281), 286 (298), 317 (324), 348 (368), 
(379.5), 412 (414), 449 (450), 457.5 (460), 539 (534), and 
548 (544) cm^^, involve both the movements of the Ca 
atoms and the MnOe octahedra. Therefore, the distinc- 
tion between external and internal modes appears to be 
less clear in CaMnOa than in PrMnOs. It follows that the 
Ca vibrations may participate significantly in the dipole 
moment of some of these modes. For instance, the dipole 
moment of the Bi„ mode at 266 (284) cm~^ (see Fig.^ is 
assumed to be mainly caused by the z component of the 
Ca vibrations and by the active 02 vibrations, namely 
the stretching and in-phase bending modes. Similarly, 
the B3„ mode at 340 (345) cm"! (see Fig.^ — in which 
all the atoms oscillate along the x direction with com- 
parable amplitudes — is expected to produce a rather 
significant dipole moment, which should be enhanced by 




FIG. 5: (Color online) Movements of atoms for a selection of 
Biu, B2u, and B^u modes of PrMnOs with Pnma structure. 
The symbols and the direction of the axes are defined in Fig. [4] 

the Ca and Mn displacements in the opposite direction 
from that of the oxygen atoms. In contrast, the Ca mo- 
tions do not play any essential role in the dipole mo- 
ment induced by the other modes at (281), 286 (298), 
361 (364), 539 (534), and 548 (544) cm^^. In particu- 
lar, in the B3„ mode at 361 (364) cm^^ (see Fig.[6| the 
modulus of calcium displacement vectors happens to be 
relatively large. However, their x component responsi- 
ble for the dipole moment is small. Finally, as shown in 
Fig.|6]the modes at 276 (275), 332 (346.5), and 510 (504) 
cm"^— as well as those at 494 (496), and 584 (613) 
cm^^ — mainly correspond to vibrations of the MnOg 
octahedra. So, these modes in which the Ca atoms can 
be considered to be at rest, can be classified as purely in- 
ternal. They consist of complex motions of the Mn and 
O atoms, except for the two modes at 494 (496), and 584 
(613) cm~^ where the displacements of the 02 atoms cor- 
respond to nearly pure stretching vibrations (~85%) , and 
pure out-of-phase bending vibrations, respectively. Nev- 
ertheless, although the amplitudes of the displacements 
of the Mn and O atoms are comparable, the dipole mo- 
ment of these modes should be mainly due to vibrations 
of the O atoms, apart for the two at 276 (275), and 332 
(346.5) cm~^ (see Fig.[6| where the Mn and 01 motions 
play similar roles. More generally, the dipole moments 
induced by the modes of frequency greater than 400 cm^^ 
are mainly due to movements of the O atoms. For the 
mode at 412 (414) cm~^ the 02 vibrations dominate. For 
those at 457.5 (460), 494 (496), 539 (534), and 548 (544) 
cm""'^ the motions of the 01 atoms play the main role. 
Only for the mode at 510 (504) cm~^ (see Fig.[6| the 
active components of the displacements of the 01 and 
02 atoms are large and comparable, but in the opposite 
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B^^ (284) 63^ (345) 63^ (364) 




FIG. 6: (Color online) Movements of atoms for a selection of 
Biu, B2u, and B^u modes of CaMnOs with Pnma structure. 
The symbols and the direction of the axes are defined in Fig.|4] 

direction. 

It follows that the vibrational patterns of IR phonon 
modes in PrMnOa and CaMnOa should be quite different 
between them, except for the modes at high frequencies. 
In fact, among the latter, only the i?3„ modes at 580 
(584) cm-i in PrMnOg (see Fig.[5| and 584 (613) cm"! in 
CaMnOs correspond to the same atomic displacements. 
Note that this phonon mode — involving mainly mixed 
vibrations of Mn and 02 atoms along the y axis (see 
Fig. [5]) — is also obtained in LaMnOa (at 506 cm"-'^) by 
means of the LDC reported by Smirnova [53 . Therefore, 
it is plausible to expect that this By,^ mode could be com- 
mon to the fJMnOa compounds with Pnma structure. 
Note also that the expected small value for the oscilla- 
tor strenght of this phonon mode is consistent with the 
experimental data. On the other hand, like the and 
B3„ modes at 571 (562) and 557 (555) cm"! in PrMnOg, 
the vibrational patterns of the and B^u modes at 548 
(544) and 539 (534) cm~^ in CaMnOa are also related to 
one another by a rotation of 90° around the y axis. As 
mentioned above, the frequency difference A/ between 
the corresponding Biu and B-^^ modes can be seen as 
a measure of the deviation of the orthorhombic struc- 
ture from the ideal (cubic) perovskite structure. There- 
fore, since the structural distorsion effects in PrMnOa 
are stronger than those in CaMnOa, we expect to ob- 
serve AfpMO > AfcMO, which is in good agreement 
with the experimental data. Finally, it should be no- 



ticed that the vibrational patterns of the Biu and B^u 
modes at 548 (544) and 539 (534) cm"! in CaMnOa in- 
volve small amplitudes for both of these modes, which is 
in poor agreement with the experimental data (see Table 
|l]). However, taking into account the complexity of struc- 
tural properties of these materials, we believe that the 
overall agreement between experiment and theory should 
be considered as reasonable. 



IV. CONCLUSION 

We have studied the IR reflectivity spectra of poly- 
crystalline orthorhombic PrMnOa and CaMnOa, the par- 
ent compounds of charge ordering manganites. LDC has 
been performed using a shell model for the two materi- 
als. The calculations have been carried out while opti- 
mizing the model parameters in order to match the cal- 
culated structural parameters, Raman, and IR phonon 
frequencies with their experimental values. By mak- 
ing a comparative study of the experimental and calcu- 
lated frequencies we have been able to find a reasonable 
assignment of the IR modes observed in PrMnOa and 
CaMnOa. In agreement with the calculations, we find 
a close relationship between the IR phonon frequencies 
of the two materials. This result shows that, contrary 
to what might be expected in an harmonic oscillator ap- 
proximation, the frequencies of medium- and high-energy 
phonons in PrMnOa are higher than in CaMnOa. Nev- 
ertheless, the absence of any correlation between the ex- 
perimental amplitudes indicates that the spectra of the 
two materials are rather different. A comparative analy- 
sis of calculated atomic displacements allows to account 
for the low-energy IR modes which appear in the spec- 
trum of CaMnOa as weak peaks. More generally, the 
present LDC shows that most of the IR modes differ sig- 
nificantly from PrMnOa to CaMnOa, and correspond to 
rather complicated atomic vibrations that cannot be as- 
signed to one type of vibration (external, bending, or 
stretching modes). So, as expected, it emerges that the 
two compounds have a different structure; nevertheless, 
CaMnOa does not seem to have a structure much closer 
to the ideal cubic perovskite structure than PrMnOa 
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